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Abstract. - A recent exponent inequality is applied to a number of dynamical growth models. 
Many of the known exponents for models such as the Kardar-Parisi-Zhang (KPZ) equation are 
shown to be consistent with the inequality. In some cases, such as the Molecular Beam Equation, 
the situation is more interesting, where the exponents saturate the inequality. As the acid test 
for the relative strength of four popular approximation schemes we apply the inequality to the 
exponents obtained for two Non Local KPZ systems. We find that all methods but one, the 
Self Consistent Expansion, violate the inequality in some regions of parameter space. To further 
demonstrate the usefulness of the inequality, we apply it to a specific model, which belongs to a 
family of models in which the inequality becomes an equality. We thus show that the inequality 
can easily yield results, which otherwise have to rely either on approximations or general beliefs. 



Growing interfaces is an important topic in Statisti- 
cal Physics since the mid 80s, booming after the seminal 
works of Edwards and Wilkinson (EW) [1] and Kardar, 
Parisi and Zhang (KPZ) 'S^. Since then, there has been 
great interest, both theoretical and experimental in such 
systems. This success is partially related to the fact that 
these systems exhibit generic properties of dynamical sys- 
tems out of equilibrium, such as dynamical steady states, 
and dynamical phase transitions 

Nevertheless, in spite of the tremendous progress, ex- 
cept for a number of one dimensional exactly soluble prob- 
lems [UEllSHl] and one two-dimensional problem |5, the 
sets of exponents given in the literature for many systems 
belonging to this class, vary considerably from author to 
author and depend strongly on the method of derivation 
[2Hill5HT5]. Interestingly, simulations suffer from a simi- 
lar problem, namely a large variety of results, which do 
not agree with each other - and even less agreement is 
found between them and the theoretical predictions [TB] . 
Researchers do not even agree as to whether an upper 
critical dimension exists or not [TT | [T4 l [T5 l [T7] . This is very 
different from the situation in equilibrium phase transi- 
tions, where methods as different as high temperature ex- 
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pansion, momentum space RG and real space RG yield 
very close exponents |18j . and those being in agreement 
with simulations as well as experimental results. Under 
such circumstances rigorous results that can put bounds 
on the exponents describing the system are obviously most 
valuable. 

Recently, a general response-correlation inequality in 
dynamical systems has been derived [19j . This inequal- 
ity is an extension of the Schwartz- SofFer inequality devel- 
oped originally for quenched random problems |20j . The 
inequality is of a generic nature and relates the response 
at steady state of some measurable physical field to an ex- 
ternal disturbance, to the time dependent correlations of 
that physical field. In the following, we explore the con- 
sequences of this inequality in various growth models. We 
begin with a short summary of the inequality, followed by 
a discussion of various analytical results that violate the 
inequality. This shows that the inequality is strong enough 
to point out problems in various approximation schemes, 
which may be the first step leading either to improving or 
discarding those schemes. Last, we show that in certain 
cases the inequality implies an equality, and demonstrate 
its potential usefulness. 

Many interesting growth models, and actually many 
other dynamical systems can be described in terms of some 
physical field, h{r, t) driven by a "noise" field, ri{r, t). This 
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very broad setup includes any of the Langevin field equa- 
tions, such as growth models of the KPZ family [5] and 
its many variants [3l|4] , noise driven Navier-Stokes pTll22] 
etc... but even more complicated cases where such an 
equation cannot be explicitly written. In the context of 
growth models, the physical field h{r, t) is simply a height 
function. Strictly speaking, h(r, t) depends not only on the 
noise field at earlier times but also on initial conditions. 
The dependence on initial conditions decays, however, in 
time, and we are left with an implicit relation between the 
Fourier transform of the field and the Fourier transform 
of the noise /i(q, w) = /i{q, oj; 77(1, cr)}, where ?7(1, cr) is a 
Gaussian random field with zero mean (77(1, a)) = 0, and 

(77(1, (7)77(01, 0) - 2i?o(^ (r)5{\ + m)5(CT + <;) . (1) 

A comprehensive description of a physical system can be 
achieved in terms of the linear response function, which 
captures the response of the system to an external pertur- 
bation, and the 2-point correlation function, which cap- 
tures the spatio-temporal structure of the physical field. 
The response function, G {q,uj), is defined by 



(2) 



and the 2-point correlation function, ^{q,uj), is defined 

by 

{h{q, Lj)h{-p, -a)) = c^)^(q - p)S{uj - a) . (3) 
In Ref. [19^ the following exact inequality was proven 



2\G(q,UjfD^{q,u:)<^q,L0). 



(4) 



This general inequality can be turned into a powerful expo- 
nent inequality. To achieve that let us recall the standard 
dynamic scaling picture in surface growth (see for exam- 
ple [3l|4l|25]) and the corresponding exponents. First, the 
equal time 2-point correlation function 



A(g) = $(<z,i = 0) 



dcLi$(q, u) , 



(5) 



describes the static properties of the growing interface. 
It generally holds that A(g) behaves as a power law in q 
for small q, namely A(g) oc . In the context of growth 
models F — d+2a where d is the substrate dimension, and 
a is the roughness exponent [3] . The larger the roughness 
exponent a, the rougher the interface is. 

The characteristic frequency, uJc{q), associated with the 
decay in time of the correlation is given by [TOilTS] 



7r$(g,t = 0)/A(q) oc q 



(6) 



where z is known as the dynamic exponent. ioc{q) has 
units of inverse time, and therefore, describes a charac- 
teristic equilibration time tx of the system as a function 



of its size L, namely tx oc i^. Therefore, the larger the 
dynamic exponent, the larger the relaxation time is, and 
the slower the dynamics. 

Another characteristic frequency, which can in general 
be different from cociq), can be obtained from the response 
function, i.e. uj^^{q) — G{q,0), which for small q behaves 
like 

^i?'(9) = G(g,0)ocg-^ (7) 

The new exponent z describes the relaxation time of the 
system in response to an external perturbation. 

Last, in this work we will focus our attention on the sys- 
tems driven by spatially (but not temporally) correlated 
noise. In other words, we will consider the family of bare 
spectral functions Do{q,ijj) (i.e., the noise correlators in 
Eq. ([1])) that for small q and uj have the form 



Do{q,u) = Bq-^" with cr > 



(8) 



The required exponent inequality is obtained now by 
setting a; = in Eq. (jU, 



2z + 2(j <r + z. 



(9) 



The equation above relates three independent expo- 
nents. It turns out, however, that most of the stochastic 
dynamical systems studied in the literature belong to one 
of two, not necessarily mutually exclusive, classes 19 . In 
each of those classes the exponent z obeys a different scal- 
ing relation, relating it to the other exponents such as 
F, z and cr. (Note that we do not claim that all stochas- 
tic dynamical systems belong to one of the two, e.g. the 
stochastic Sine-Gordon systems p3l[24] ). 

Class I is that of generalized Hamiltonian systems where 
S = F — 2(7. This is the case for example in Model A of Ho- 
henberg and Halperin [18] (i.e., the celebrated (j)'^ theory) 
where z = F. Class II is the class of Galilean invariant sys- 
tems. It is composed of systems whose dynamics is given 
by a Langevin field equation of the form 



(10) 



supplemented with the condition SFq{h} /Sh (q = 0) = 0. 
For growth models of the KPZ family [3 this corresponds 
to the fact that the dynamics is invariant to translations 
in the h-direction, i.e. with respect to h ^ /i -I- Const, also 
known as the Galilean Invariance. For this class it was 
shown in Ref. |19j that z = z. It is important to mention 
that the proof is based on the Family- Vicsek dynamical 
scaling [31[ni[T3[Tli25j 



q 



(11) 



Note, that most of the growth models belong to Class 
II, including the KPZ equation [2], the Molecular Beam 
Equation (MBE) [26^ and many more [3j|4]. Actually 
the noise driven Navier-Stokes equation [3T1[52], various 
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wetting-front models [27H1S] and crack propagation equa- 
tions [301131] also belong to this class. 

Using the relation z — the general inequality ^ re- 
duces to the simpler 



< r - 2ct < r , 



The above inequality is correct for any set of exponents, 
whether strong coupling or weak coupling. In fact, for 
weak coupling exponents, we expect the leftmost inequal- 
ity to hold as an equality. 

The most famous systems belonging to Class II are 
growth models of the KPZ [1|3] and MBE families 
[3l[26l|32l[33] . In those systems, we can obtain more than 
just an inequality relating the dynamic exponent z and 
the steady state exponent F. This is because the two are 
related by a scaling relation (F — (i)/2 + z = 2 for the 
KPZ family and (F - d)/2 + z = 4 for the MBE family 
[3l[26]. (Recall that a = (F — (i)/2 is the roughness expo- 
nent). It is easy to verify that all the known results for 
regular KPZ obey the inequality, which becomes in this 
case 2 < (d -t- 4)/3. For d = 1 all the analytic methods 
recover the exact results for regular KPZ, namely z = 3/2 
that obeys the inequality z < 5/3. Therefore, there is no 
surprise that in one dimension the inequality is obeyed for 
regular KPZ by all the methods. For higher dimensions, 
although the exponents obtained analytically 
deviate most considerably from the simulations |16| all 
methods yield T > d and z < 2, and so the inequality 
becomes less tight. 

For the MBE equation (5^ , the inequality becomes 
z < (d-|-8)/3. The one-loop Dynamic Renormalization 
Group (DRG) result [26|, as well as the Self Consistent 
Expansion (SCE) [34] yield z = {d + 8)/3, and so the 
inequality is saturated. 

To demonstrate the usefulness of the exponent inequal- 
ity, we use it to test various theoretical approaches to the 
study of the KPZ family, and concentrate on two non-local 
extensions to the KPZ equation. From the early days of 
the field it was clear that the simple models such as the 
EW and KPZ models show systematic deviations from ex- 
perimental data. However, since the general scaling pic- 
ture developed in this context [3]l4p5] usually applies, var- 
ious models with correlated noise f55^ conserved [2Sll32l33j , 
nonlocal [37H39] and fractional dynamics [7l|40] have been 
developed to account for important physical effects that 
were not taken into account beforehand. 

A Non-local KPZ (NKPZ) equation has been introduced 
in [37] to account for the non-local hydrodynamic interac- 
tions in deposition of colloidal particles in a fluid. It was 
later generalized to spatially correlated noise {a ^ 0) in 
Ref. [35]. The non-local KPZ equation suggested in [37] 
for the height function, h{r, t) of the deposited material is 
given by 



dh (r,t) 
dt 



(13) 



+ J dr'g (r') \/h (r + r', t) • V/i (r - r', t) + ?] (r, t) , 
where the kernel g (r) has a short range part, XqS (r), and a 



(12) long-range part ^ Xpr^ . In Fourier space, g (q) = Aq 



Xpq^'' {p is the non- locality parameter). For simplicity, 
we discuss here only the case Aq = 0. The noise has zero 
mean, but is allowed to have general correlations of the 
form given by equation (jS]). The strong coupling solution 
found by the DRG is 



zdrg 



(d - 2 - 2p) (d - 2 - 3p) 
{3 + 2-P)d-6-9p 



(14) 



Because of the extra scaling relation between a and z [371 
138] , namely a -I- z 2 — p, there is only one independent 
exponent in NKPZ. 

The above result violates the inequality (IT^ over a 
whole range of parameters defined by Fd j^g — zdi^q^'^^^ < 
0. To be more concrete, let us specialize to the case a — 
where the inequality, TjyfiQ — zd^g < is violated in 



d~2 



< p < Po (d) for d < do . 



and 



po{d) < p< 



d~2 



for d > do ■ 



(15) 



(16) 



3.395 is the 



(2-d 
i2T~dln2 I 

solution of po(do) = and W (x) is the Lambert func- 
tion. The shaded region in Fig. [T] is the region in the (d, p) 
plane where the inequality (i.e. Eq. with cr = 0) is 
violated. 
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Fig. 1: Violation of the response-correlation inequality (|12[l 
with (J = by the DRG method derived in Refs. [3711381 occurs 
in the shaded area enclosed by the curves p = (d — 2)/2 and 
p = po (d) in the phase diagram. 

We now turn to results derived for the same NKPZ 
model (fT4)l using an improved Mode-Coupling (MC) ap- 
proach originally derived for the local KPZ problem in |14) 
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and later applied by Hu and Tang [IT] to the non-local 
case. In Fig. [5] we give the solution for the dynamic ex- 
ponent z for d = 1,2,3 reported in [41]. The shaded area 
bounded by the line z = (d -I- 4 — 2p) /3 marks the region 
where the inequality is violated by the MC prediction. 
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Fig. 2: (Colour online) The Dynamic exponent z as a function 
of the non-locality parameter p using the mode coupling result 
(solid (blue) line - reproduced from Ref. [41]). The part of 
the solid (blue) line within the shaded region is excluded by 
the inequality p2|) with ct = 0. As can be seen the mode 
coupling result violates the inequality in certain regions in all 
dimensions. 



Next, we check for violations of the inequality in results 
obtained by Tang and Ma [32] generalizing the Flory-type 
Scaling Approach (SA) of Hentschel and Family [9| to the 
non-local case. The strong-coupling dynamical exponent 




The nonlocality parameter p 



Fig. 3: Violation of the response-correlation inequality p2|l by 
the SA method derived in Ref. 42] occurs in the shaded area 
enclosed by the curves p = — 1 and a = {d + 1 + p)/2 in the 
phase diagram. 



obtained by using this method is 

(2-p) {d + 2-2(j) 



ZSA 



d + 3-2(T 



(17) 



It turns out that this solution violates the inequality in a 
whole region of parameter space defined by 



d d+l 

:t < o" < TT 

2 2 



and 



1 



p d 
-<"<2 



for p > — 1 



for p < — 1 



(18) 



(19) 



These results are presented graphically in Fig. [S] 

Last we tested the results for NKPZ of the Schwartz- 
Edwards Self-Consistent Expansion (SCE) [THKIS. This 
method predicts a whole zoo of possible phases all of which 
are consistent with the inequality, and therefore the only 
method so far that does not contradict the inequality. 

To check if the picture we get so far may be more gen- 
eral we consider another non local system, which has been 
introduced and analyzed in Ref. [35]. In this model, the 

nonlinearity in Eq. is replaced by J dr' '^''^jrLJ'^'i-p ^ i 
such that the contribution to the growth at each point 
comes from the interaction of the gradient at that point 
with all the other points on the interface, rather than 
from all pairs at equal distance. The Scaling- Approach 
mentioned above yields the same results as before for this 
model as it depends only on the dimensions of the various 
quantities involved, and ignores the exact spatial struc- 
ture, and therefore, gives rise to violations of the inequal- 
ity here too. The situation with the DRG approach is even 
more severe as the theory generates more relevant terms 
under renormalization, and therefore not surprisingly pro- 
duces inconsistencies with the inequality. Interestingly, 
the results of the SCE method (reported in Ref. [3S]) are 
not only consistent with the exact one-dimensional result 
[5] but the exponents in all the phases are consistent with 
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-1 
The non-locality parameter p 




-1 
The non-locality parameter p 




-1 
The non-locality parameter p 

Fig. 4: (Colour online) The dynamic exponent z as a function 
of the non-locality parameter p for uncorrelated noise {a = 0) 
in d = 1,2,3 dimensions resulting from the Self-Consistent 
Expansion f3Sl . For d = 1 the solution saturates the inequality 
for —1/2 < p < — 1/5 . Note that for d — 2,3 two phases 
are possible in some range of the parameters as seen from the 
fact that 2 is a multivalued function of p (solid (blue) line 
weak coupling solution, and dashed (red) line strong coupling 
solution). Also note that for d = 2 and p = there is only one 
phase as the solid (blue) line does not exist for p = 0. 



our stronger, cr-dependent inequality over the whole pa- 
rameter space and all dimensionalities. In Fig. |3]we give 
the SCE solution for the dynamic exponent z for d = 1, 2, 3 
for uncorrelated noise {a = 0)(but as stated above the in- 
equality holds also for a > 0). 

The inequality is shown to be a useful tool in detect- 
ing shortcomings of various analytical methods applied to 



the NKPZ family. This suggests that results obtained by 
those methods should be suspected even in cases where 
the inequality is not violated ,at least until the origin of 
violation of the inequality is understood. 

The last demonstration of the usefulness of the inequal- 
ity, is when a system belongs to both Class I and Class II, 
i.e. when the system is Hamiltonian as well as Galilean 
Invariant. In these cases the inequality becomes an equal- 
ity [IH], namely z = T — 2a. This reduces the number of 
unknown exponents by one, and can be a powerful tool 
when dealing with certain nonlinear growth models. For 
example, a surface evolving due to deposition in the pres- 
ence of surface tension can be described by the Hamilto- 
nian H ^ ^ 



Jd'^rJl + (V/i) , which leads to the follow- 



mg equation 



Vh 



'i + ivhY 



(20) 



This model is believed to belong to the EW universality 
class [3] on the basis of symmetry arguments, or by us- 
ing a small-gradient expansion that leads directly to the 
EW equation [1]. The inequality shows rigorously that 
z — d + 2a, avoiding issues of the validity of the small- 
gradient expansion. This actually allows deriving quick 
and accurate results when applicable. 

To summarize, in this paper we have shown how to use 
a recent inequality derived in [19] to growth models de- 
scribed in terms of stochastic field theories. We show that 
the inequality, which involves the correlation and the re- 
sponse functions can be translated into a simple inequality 
for the scaling exponents F, z and the noise correlation ex- 
ponent a, in cases where a single q dependent time scale 
exists. Although being extremely simple, this inequality 
can be quite powerful when examining analytical, numer- 
ical and experimental results. 

To demonstrate the utility of the inequality, we reviewed 
analytical results for two non-local KPZ models, obtained 
by using four different methods: Dynamical Renormal- 
ization Group, Mode-Coupling, Scaling- Approach and the 
Self-Consistent Expansion. Interestingly, the first three 
methods yield results which contradict the inequality for 
a whole set of parameters, while the Self- Consistent Ex- 
pansion is the only one which never violates it. This has 
an important implication on the choice of analytical tools 
when dealing with such stochastic models. Also, we have 
pointed out the observation that whenever the system is 
Hamiltonian as well as Galilean invariant, the general in- 
equality implies an equality z = T — 2a, and can be quite 
powerful when applicable. Last, even in cases where the 
inequality is not saturated nor implies an exact equality, 
in effect it can provide a reasonable estimate for the scal- 
ing exponents. In the cases we have reviewed, the bounds 
implied by the inequality do not deviate a lot from the real 
result, especially in low dimensions. This suggests using 
the inequality to derive quick approximate estimates for 
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nonlinear systems. 

This work has interesting experimental implications. If 
in a given experimental system z 7^ z we get a strong 
constraint on the kind of theory that could describe the 
phenomenon under consideration, namely a theory that 
cannot have Galilean invariance. Similarly, if experimen- 
tally r ^ z the theory cannot be Hamiltonian. This can 
provide useful guidelines for both theorists and experimen- 
talists who may wish to model the system. 

There are obviously some interesting open questions. 
For example, are there relevant growth models that vio- 
late the Galilean invariance, and therefore have distinct z 
and z? More generally, the question here is regarding the 
existence of interesting classes other then the two identi- 
fied so far. Another direction could be to explore similar 
inequalities for systems described by more than one-field 
[43] . and systems with quenched disorder [44] . 
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